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ABSTRACT 



The communication channel between a satellite and a ship earth station (SES) is 
described by a model which includes multipath fading, doppler shift and noise. Multi- 
path fading is caused by reflections from the sea surface. These reflections can affect the 
system performance, especially at low elevation angles or when SES is using low gain 
antennas. Doppler shift is a very important effect when using low altitude satellites, be- 
cause of the high velocities involved. 

This thesis describes and presents a software simulator for multipath fading in the 
maritime communications environment. Analysis of throughput of an unslotted Aloha 
maritime mobile satellite communication channel is also presented. 
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I. INTRODUCTION 



Future satellite systems for maritime, aeronautical and land mobile applications 
(with mobile terminals for low data rate and digital voice transmission) will have small 
antennas having gains between 2 and 12 dB and a broad beamwidth. These antennas 
pick up multipath reflections of satellite signals causing severe fading of the received 
signals. 

Several fading simulators, using hardware, software, or both, have been designed. 
Among them are: JPL, CRC and Virginia Tech. All of them deal with land mobile 
communications and with geostationary satellites. 

In the future, low altitude satellites will become important alternatives for mobile 
communications because of the crowded geostationary orbit. Studies have demonstrated 
that five satellites in circular orbit will provide continuous single coverage of the earth; 
seven satellites will provide double coverage at some minimum elevation angle. By the 
same criterion it was found that eleven satellites can provide continuous triple coverage 
and fourteen satellites, quadruple coverage. Continuous single coverage means that any 
two stations on earth can communicate using this constellation of five satellites at any 
moment. Double coverage, on the other hand, will provide an additional communication 
path; triple coverage will provide three diflerent paths and so on. 

Besides the economical achievement of maximum coverage, many other factors have 
to be taken into account in system design. These results may prove of some interest to 
planners of satellite systems for mobile communications and surveillance [Ref. 1] 

Low earth orbits (LEO) are generally circular or slightly eccentric. In order to 
simplify the analysis and equations described in this thesis, the focus will be on circular 
orbits. 

This thesis presents a fading simulation program for mobile maritime satellite com- 
munications involving low altitude satellites. 

Chapter II describes the different signals involved in the maritime communications 
channel, the atmospheric effects that affect them, and also demonstrates the application 
of the theory of Rician fading statistics. 

Chapter III describes the simulation program, explains the purpose of each of its 
components and illustrates some results obtained by running several simulations for 
different values of frequency and satellite's velocity. 
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Chapter IV describes the throughput of unslotted Aloha with power capture prob- 
ability in a Rician fading environment. 

Finally, Appendix A contains a listing of the computer program of the simulator, 
and Appendix B shows the derivation of the formulas used in chapter IV. 
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II. PHYSICS OF MARITIME MOBILE SATELLITE PROPAGATION 



A. INTRODUCTION 

Any signal in maritime mobile communications is typically dominated by a multi- 
path fading process caused by reflections coming from the surface of the sea. These 
multipath signals are expected to be stronger at low satellite elevation angles, where the 
main lobe of the antenna of the ship earth station (SES) picks them up. No specular 
component is expected, except for the very calm, mirror-like sea; therefore, a model with 
Rician statistics is suggested in this chapter. 

B. MULTIPATH FADING AMPLITUDE 

It is well known that multipath reflections from the sea surface consist of specular 
and diffuse components. The relative strength of these components is due to wave height 
and to elevation angle, with the former the dominant factor [Ref.2]. For a smooth sea 
state the specular component dominates. As wave height increases, the diffuse compo- 
nent also increases, and over a wave height of 0.5 m, the diffuse component becomes 
constant. Sea state probability studies have been done in the past and it was concluded 
that wave heights are usually in a 0.5 to 5 m. range [Ref. 3]. Based on these observations, 
it is assumed that sea conditions will be rough most of the time and thus amplitude 
fading can be expressed in statistical terms. 



C. STATISTICAL MODEL FOR MULTIPATH FADING AMPLITUDE 

For a satellite maritime link, the dominant multipath component is the diffuse 
component. The signal can therefore be represented by two components: direct and dif- 
fuse (see Figure 1). 
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The direct component is the line of sight (LOS) signal from the satellite. The diffuse 
component is the resultant of many random scattered waves arriving in a uniform an- 
gular distribution (Ref. 4]. If each of the waves is represented by a random phasor, the 
total diffuse component at the output of the antenna is 



R dif = r expfrtf] = expC/0/] (2.1) 

/=i 



where 

• r is the amplitude of the diffuse component 

• 0 is the phase of the diffuse component with respect to the direct component 
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• A, is the amplitude of the i"' scattered wave with respect to the direct component. 

If the scattered waves are sufficiently random and their phase is uniformly distrib- 
uted over 2n, Beckmann shows that the diffuse component will be Rayleigh distributed 
with a density function given by [Ref. 5] 

pir) = exp[^ ~ V ] r ^ 0 ( 2 - 2 ) 



p(r ) = 0 r < 0 



where 

• r is the received diffuse signal level (voltage) 

• a = E{r 2 } 

The Rayleigh distribution requires the scattered waves to have a uniform phase dis- 
tribution over an interval of 2 n. If the terrain does not provide for this requirement, the 
Rayleigh distribution will not be valid [Ref. 6]. 

The vector sum of a direct signal, with a Rayleigh fading component, results in a 
composite signal with Rician envelope statistics [Ref. 7], and the probability density 
function is: 



p{r) = 



ir 

a 



exp - 



'- 2 + A\ 
a 



Ko( 



2A 0 r 



r > 0 



• (2-3) 



p{r) = 0 



r < 0 



where 

• r represents the envelope of the Rayleigh fading signal 

• a = E{r 2 } 

• A 0 is the amplitude of the direct signal 

• I 0 is the modified Bessel function of zero order 

By normalizing the direct component to unity, this reduces to: 



p(R) = 



2 R 

k 




R 2 + 1 
k 




R> 0 



(2.4) 
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P {R) = 0 



R< 0 



where 



• k is equal to the power in the multipath component over the power in the direct 
component, that is: 




• R is the ratio of the amplitude of the fading signal to the amplitude of the direct 
component 




The parameter k is called the Rice factor. It is considered a key parameter in the 
Rician distribution, and its decibel equivalent is K = 10 log 10 A. Figure 2 illustrates se- 
veral Rician distributions for different values of A 0 . 

The phase of the received signal is not uniformly distributed [Ref. 5] and its density 
function is given by: 



P(°) = — ex P 




(1 + G v 7r exp(G 2 )(l + erf[G))) 0<d <2 n (2.5) 



where 



• G is defined as: 



G = 



A 0 cos 6 



erf (G) is the error function defined as: 
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Figure 2. Rician density function for a = 1, f = 869 MHz, V = 0 Kts. 

1. Cumulative Distribution Function (CDF) 

The cumulative distribution function (CDF) refers either to the function G(A) 
for which the signal exceeds a level A, or to the function F(A) for which the signal is 
below a level A. The former is the one used in this thesis. 

The theoretical expression for the Rician CDF is obtained by integrating 
equation (2.3): 



G(A) = P{r>A) = \ p{r)dr (2.6) 

Figure 3 illustrates a Rician CDF. 
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RICIAN CDF 




Figure 3. Cumulative Distribution Function (CDF) 

Usually the CDF is plotted in Rayleigh paper which has the especial character- 
istic that a Rayleigh curve drawn in it is a straight line. This is useful for evaluating any 
outcome by visually comparing it with the Rayleigh line. Figure 4 illustrates the same 
Rician CDF seen in Figure 3 but on probability paper. 
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RICIAN CDF 




-15 -10 -5 0 5 

RECEIVED SIGNAL RELATIVE TO LINE OF SIGHT (dB) 



Figure 4. Cumulative Distribution Function (CDF) on probability paper. 

D. ATMOSPHERIC EFFECTS ON THE FADING SIGNAL COMPONENTS 

The direct and the difTuse components are both affected by Ionospheric and 
Tropospheric efTects. These effects include Faraday rotation, group delay, absorption, 
dispersion, refraction and scintillation, all of which result from interaction with the 
earth's magnetic field and the ambient electron content as the signal passes through the 
Ionosphere. Table 1 [Ref. 8] summarizes some typical maximun values obtained for an 
elevation angle of 30 °, one way propagation, at two different frequencies. Faraday ro- 
tation can be minimized using circular polarization. Scintillation can be important at low 
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latitudes (9° N to 21 0 S) and auroral latitudes (above 59 ° X), but they are insignificant 
for elevation angles greather than 10 ° and frequencies below 10 Ghz [Ref. 9]. 



Table 1. ESTIMATED MAXIMUM IONOSPHERIC EFFECTS [REF. 8J: For 

an elevation angle of 30 ° , one way propagation, and a zenith electron 
column of 10 18 electrons/m J 



Effect 


Frequency 

dependence 


Magnitude 


850 MHz 


1600 MHz 


Faraday Rotation 


1 if 


o 

o 


42 0 


Group delay 


1 if 


0.35 s 


0.1 s 


Variation in direction of arrival 


1 if 


16 sec of arc 


4.7 sec of arc 


Refraction 


1 if 


> 50" 


> 14" 


Absorption (mid-lats) 


1 if 


> 0.014 dB 


> 0.004 dB 


Dispersion 


1 'f 


0.65 nsec/ MHz. 


0.1 nsec MHz 



Tropospheric effects result from moisture in the atmosphere and tend to become 
more severe as frequency increases. Table 2 [Ref. 10] summarizes some estimated values 
for an elevation angle of 30 °, one way propagation, at two different frequencies. 

The antenna pattern will influence the magnitude of the diffuse component. The 
antenna gain rolls ofT below the horizon; therefore, the contribution of the diffuse com- 
ponent is due to scattered waves arriving from angles above the horizon. 

The current treatment assumes that the diffuse component is 8 to 14 dB (mis) below 
the direct component [Ref. 9]. 
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Table 2. ESTIMATED TROPOSPHERIC ATTENUATION [REF. 10]: For an 

elevation angle of 30 °, and one way propagation 



Effect 


Magnitude (dB) 


850 MHz 


1600 MHz 


Clear air absorption 






3 g'm 3 (dry) 


0.06 


0.07 


7.5 g/m 3 (average) 


0.06 


0.07 


17 g'm 3 (moist) 


0.06 


0.07 


Cloud attenuation 






0.5 g'm 3 , 1 Km thick 


< 0.01 


< 0.01 


1 g/m 3 . 2 Km thick 


< 0.01 


< 0.01 


Fog attenuation 






0.05 g m 3 (average). 0 to 75 Kmt ht. 






0.05 g'm 3 (heavy). 0 to 150 mt ht. 






Rain attenuation 






5 mm h 


< 0.01 


< 0.01 


25 mm h 


< 0.1 


< 0.1 
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III. SOFTWARE FADE SIMULATOR FOR A MARITIME MOBILE 

COMMUNICATION CHANNEL 



A. INTRODUCTION 

Arredondo et al. [Ref. II] describe a hardware simulator to simulate multipath fad- 
ing for mobile radio, and Karim [Ref 12] describes a different solution to implement this 
simulator using software. This thesis follows both reference models, adding a constant 
signal (direct component) to obtain Rician statistics. It emphasises the doppler effect 
caused by the high velocity of low altitude satellites and the vehicle's velocity. Figure 5 
is a block diagram of the suggested model. 



B. SIMULATOR DESCRIPTION 

1. Rayleigh Multipath Fading Generator. 

The most important part of the simulator is the Rayleigh multipath fading com- 
ponent. Its implementation follows the model of Arredondo et al. [Ref 11] and Karim 
[Ref. 12]. Two independent Gaussian sources with zero mean and equal variance are 
generated. They are passed through a shaping filter and finally added in quadrature to 
obtain a signal with an envelope that is Rayleigh distributed in amplitude and uniformly 
distributed in phase. 

The generator is illustrated in Figure 6. A key aspect in the design of this gen- 
erator is the shaping filter. 



12 



n Ricinn fading model 




RICIAN 

FADING 

SIGNAL 



Figure 5. Rician fading simulator block diagram 
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A RAYLEIGH FADING SIMULATOR 




FADING 

SIGNAL 



Figure 6. Rayleigh fading generator block diagram 



Shaping Filter: The spectral density S(a>) of the envelope of the Rayleigh fading 
signal received by an omnidirectional antenna is given by [Ref. 12] : 

) 2 r°' S |o> - to* :£ to, (3.1) 

where 
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• co 0 is the carrier frequency. 

• co rf is the doppler shift due to the velocity of the vehicle and the velocity of the 
satellite. 

• a 2 is the variance of the two independent Gaussian sequences. 

Refering to equation (3.1), taking co relative to and noting that we are only 
interested in relative values of the signal, the transfer function H(co) that has to be ob- 
tained to shape the input is given by: 



H(co) = [1 - ( ) 2 ] -0 ' 25 (3.2) 

This function is unbounded for co = co d , but following the model of Arredondo 
et al. [Ref. 11], the actual transfer function has a maximum at 0.95(0,, and falls to 0 with 
a slope of -18 dB / octave. It is done this way so that the output envelope of the simu- 
lator closely fits a Rayleigh distribution. 

This transfer function will be implemented as the transfer function H(z) of a 
non-recursive digital filter with sampling frequency co, = 2co d , and the procedure is dis- 
cussed in Karim model [Ref. 12]. 



20 

//(z) = 0 z-‘ (3.3) 

1=0 

where 

• IV, is the Hamming window. 

• a, is the Fourier coefficient of the Fourier series expansion of equation (3.2). 

• z is exp(/'co7) where T = -py- . 

Figure 7 illustrates both the theoretical and actual transfer functions. 
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I 





Figure 7. a) Spectral density of complex envelope, b) Shaping filter freq. 
response.: [Ref. 11] 



Doppler Shift [Ref. 13] 

Doppler shift is caused by a length change with time in the propagation path. 
Let the signal from the satellite be A 0 cos{2nf 0 t). This signal arrives to the receiver after 
a time delay g(r), where g is a function describing the time variation in path length. 
Therefore the received signal is given by A^cos[2nf^t±g{t))]. Differentiating the phase 
we can obtain the frequency of the received wave, that is: 

-^•■|-[2^o('±*«))] =/o('±-^r) o- 4 ) 

Thus the shift frequency is: 
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(3.5) 



A/ 









D{t) 

Since g{t) = ^ where 

D{t ) is the slant range, and 
C is the speed of light, 



tfo 



di 



we have 



dg{t) d[D(t)IQ 



v 

C 



di dt 

where v is the relative velocity of the receiver with respect to the source. Due to the 
higher velocity of the satellite compared with the velocity of the receiver on earth, the 
latter becomes negligible. For the purpose of this thesis, only the radial component of 
the velocity of the satellite v, will be considered. 

Therefore equation (3.5) becomes: 



A/= i/o —q 



(3.5 a) 



For satellites 



40 = [(// + R e f + R 2 - 2 R e (H + R e ) cos 9] 



(3.6) 



where 

• H is the altitude of the satellite orbit 

• R r is the radius of earth = 6378.155 Km 

• 8 is a linear function of time for circular orbits. 

Due to the fact that earth terminals are not located at the center of the earth, 
there is a maximum Doppler shift at maximum slant range (satellite "rise" and "set"), as 
shown in Figure 8. From the same figure the following are obtained: 



[ s = R e + H 


(3.7) 


”"(t) 


(3.S) 


* r = v sin y 


(3.9) 
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where 



• R s is the total radius of the satellite orbit. 

• y is the angle described in Figure 8 

• v, is the radial velocity of the satellite, from which the maximum value of Doppler 
shift is obtained. 




Figure 8. Doppler shift from a low altitude satellite 



Therefore, the maximum value of Doppler shift is: 



v r,ma.x \ v r,max 



U max = ±fo{ J ^-) = ± 






(3.10) 



where 

• f d is the dopplcr frequency 
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• / is the wavelength of the carrier frequency 

• f 0 is the transmitted carrier frequency. 



For a circular orbit satellite, the velocity is given by: 



v 




(3.11) 



where 



He~ + - g(^ear,l) - 3.986013x10 s 



Km' 

sec 2 



From equations (3.9), (3.10) and (3.11) we obtain: 

vR 



H e 

.... N VX Re v 13.42/o 

fd,max = ±/o( ) = ±fo( FF ) - ± 



C/?, 



3 



(3.12) 



2. Direct Component part. 

Let the direct component of the signal transmitted by the satellite be 
.4 0 cos(ty 0 r). The signal at the receiver will be affected by the velocity of the satellite, the 
velocity of the vehicle, or both, causing a Doppler shift. Davarian [Ref. 14] suggests a 
hardware solution to take care of the Doppler shift in the direct component; this is 
shown in Figure 10. As a result, the signal at the receiver will be ,4 0 cos[(a> 0 + coj/] at 
satellite rise and ,4 o c0s[(<o o — t yjr] at satellite set. One of these signals is added to the 
Rayleigh fading signal to obtain the Rician fading Signal. 
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DOPPLER SHIFT 






COS(d 


>oO 

f 

— 




COS(d> d t ) 


— *6 




► COS(<o 0 + <o d )t 




MIXER 





Figure 9. Block diagram for frequency shifting of the direct path [Ref. 14 ]: co 0 is 

the angular frequency of the carrier and co rf denotes the Doppler shift. 



C. COMPUTER PROGRAM 

A computer program was developed to simulate Rician fading statistics for a mari- 
time mobile satellite communications channel. The computer program is written in 
Fortran and is included as Appendix A. 
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D. ANALYSIS OF RESULTS 



The simulation of several sequences of multipath Rician fading was done assuming 
a communication mobile maritime link with a satellite in circular orbit at an altitude of 
500 nautical miles. As was stated earlier in this chapter, the speed of the maritime vehicle 
is considered negligible in comparison to the speed of the satellite. As a result, the vehi- 
cle's speed is not taken into account. Figures 10 through 12 show the generated Rician 
sequences for values of K (Rician factor) equal to 0, 5 and 10 dB, a frequency of 869 
MHz and a velocity of 30 knots. Figures 13 through 15 illustrate the histogram of the 
preceeding sequences. Figures 16 through 31 show the cumulative distribution function 
(CDF), using probability paper, of several multipath sequences for the different combi- 
nations of the following parameters K = -26, -10, 0, 10 dB, frequency = 869, 1501 
MHz, and velocity = 0, 14357 Kts. 

In order to analize the generated Rician sequences we look at two different CDF 
graphs. Comparing Figures 16 and 28 using the same frequency and velocity for values 
of K = -26 dB and K = 10 dB, it can be noticed that for the former value of K the re- 
ceived signal is greater than the abscissa for a longer period of time. This means that for 
a larger Rice factor (K) we have more fading and therefore less strength in the received 
signal at the antenna. 

Additionally, the difference in the CDF can be noticed between the curves with 0 
velocity (0 doppler) and the curves with a velocity of 14357 Kts (max. doppler); as was 
expected, the difference is larger at higher frequencies. 



21 



K = 0 dB, f = 869 MHz, V r = 30 Kts 
FADING SEQUENCE 




Figure 10. Rician fading sequence for K = 0 dB, f = 869 MHz, V = 30 Kts. 
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K = 5 dB, f = 869 MHz, V r = 30 Kts 
FADING SEQUENCE 
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K = 10 dB, f = 869 MHz, V r = 30 Kts 

FADING SEQUENCE 
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Figure 12. Rician fading sequence for K = 10 dB, f = 869 MHz, V = 30 Kts. 
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K = 0 dB, f = 869 MHz, V r = 30 Kts 
HISTOGRAM 




R(dB) 



Figure 13. Histogram of a Rician fading sequence for K = 0 dB, f = 869 MHz, V 
= 30 Kts. 
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K = 5 dB, f = 869 MHz, V r = 30 Kts 
HISTOGRAM 




R(dB) 



Figure 14. Histogram of a Rician fading sequence for K = 5 dB, f = 869 MHz, V 

= 30 Kts. 
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K = 10 dB, f = 869 MHz, V r = 30 Kts 

HISTOGRAM 




R(dB) 



Figure 15. Histogram of a Rician fading sequence for K = 10 dB, f = 869 MHz, 

V = 30 Kts. 
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K = -26 dB, f = 869 MHz, V r = 14357 Kts 

CDF 




RECEIVED SIGNAL RELATIVE TO UNE OF SIGHT (dB) 



Figure 17. CDF of a Rician fading sequence for K 
14357 Kts. 



-26 dB, f = 869 MHz, V = 
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K = -26 dB, f = 1501 MHz, V r = 0 Kts 

CDF 




RECEIVED SIGNAL RELATIVE TO LINE OF SIGHT (dB) 



Figure 18. CDF of a Rician fading sequence for K = -26 dB, f = 1501 MHz, V 
= 0 Kts. 
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14357 Kts 



K = -26 dB, f = 1501 MHz, V r = 

CDF 




RECEIVED SIGNAL RELATIVE TO LINE OF SIGHT (dB) 



Figure 19. CDF of a Rician fading sequence for K = -26 dB, f = 1501 MHz, V 
= 14357 Kts. 
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K = -10 dB, f = 869 MHz, V r = 0 Kts 

CDF 




RECEIVED SIGNAL RELATIVE TO LINE OF SIGHT (dB) 



Figure 20. CDF of a Rician fading sequence for K = -10 dB, f = 869 MHz, V = 
0 Kts. 
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K = -10 dB, f = 869 MHz, V r = 14357 Kts 

CDF 




RECEIVED SIGNAL RELATIVE TO LINE OF SIGHT (dB) 



Figure 21. CDF of a Rician fading sequence for K = -10 dB, f = 869 MHz, V = 
14357 Kts. 
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Figure 22. CDF of a Rician fading sequence for K = -10 dB, f = 1501 MHz, V 
= 0 Kts. 
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K = -10 dB, f = 1501 MHz, V r = U357 Kts 

CDF 




RECEIVED SIGNAL RELATIVE TO LINE OF SIGHT (dB) 



Figure 23. CDF of a Rician fading sequence for K = -10 dB, f = 1501 MHz, V 
= 14357 Kts. 
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K = 0 dB, f = 869 MHz, V r = 0 Kts 

CDF 




RECEIVED SIGNAL RELATIVE TO LINE OF SIGHT (dB) 



Figure 24. CDF of a Rician fading sequence for K = 0 dB, f = 869 MHz, V = 0 
Kts. 
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K = 0 dB, f = 869 MHz, V r = 14357 Kts 

CDF 
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Figure 25. CDF of a Rician fading sequence for K = 0 dB, f = 869 MHz, V 
14357 Kts. 
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K = 0 dB, f = 1501 MHz, V r = 14357 Kts 

CDF 




RECEIVED SIGNAL RELATIVE TO LINE OF SIGHT (dB) 



Figure 27. CDF of a Rician fading sequence for K = 0 dB, f = 1501 MHz, V = 
14357 Kts. 
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K = 10 dB, f = 869 MHz, V r = 0 Kts 

CDF 




Figure 28. CDF of a Rician fading sequence for K = 10 dB, f = 869 MHz, V = 
0 Kts. 
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Figure 30. CDF of a Rician fading sequence for K = 10 dB, f = 1501 MHz, V = 
0 Kts. 
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K = 10 dB, f = 1501 MHz, V r = 14357 Kts 

CDF 




Figure 31. CDF of a Rician fading sequence for K = 10 dB, f = 1501 MHz, V = 
14357 Kts. 
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IV. CAPACITY OF UNSLOTTED ALOHA MARITIME MOBILE 

SATELLITE SYSTEMS 



A. INTRODUCTION 

Aloha is a multiple access method commonly used in data communications. This 
chapter analyzes the throughput of unslotted Aloha with power capture in a Rician 
fading environment. 

B. UNSLOTTED ALOHA THROUGHPUT 

In the Aloha system, a satellite channel with a capacity of R bits per second is 
shared by a large population of users. Each user's station randomly transmits packets 
at the channel bit rate R whenever its buffer contains one packet. When packets from 
different stations overlap (packet collision), the transmission error can be detected at the 
transmitting stations on the downlink. The stations then retransmit the packets until 
they are free from overlap. To avoid repeated collisions, the interval of packet trans- 
mission is randomized for each station [Ref. 15]. 

In unslotted Aloha there are two types of capture effect: created and natural. Cre- 
ated capture effects occur when random power levels are used. In this case, the packet 
with the higher power level may capture the receiver when its signal-to-interference ratio 
exceeds a preset threshold. Natural capture effects occur in mobile radio communi- 
cations when varying distances between mobile nodes and base stations, as well as 
channel fading, cause varying signal powers at the base station's receiver [Ref. 16 ]. 

As stated in the paper by Borchardt et al. [Ref. 17] about unslotted Aloha, let us 
assume an infinite number of uncoordinated users communicating to a base station us- 
ing fixed length packets. The channel traffic is approximated by a Poisson proccess with 
parameter G packets per packet length. Thus the probability that a tagged packet is 
overlapped by n other packets during its transmission is the arrival probability of n 
packets during the interval (r 0 — T, /„ + T) , where i 0 is the transmission instant of the 
tagged packet and T is the packet transmission time. This arrival probability is 

(2G) n 

Pr{n) = ; — exp( — 2G) (4.1) 

n\ 
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According to Borchardt et al. (Ref. 17], consider a receiver that captures a tagged 
packet in the presence of n interfering packets, if the tagged packet power sufficiently 
exceeds the joint interference power. These interferers form a stochastic proccess which 
consists of early interferers, whose transmission began in the interval (r 0 — 1, r 0 ), and late 
interferers which enter the system in the interval (/ 0 /o+ O* With n interferers there are 
2" realizations of early and late interferers. C (n) is the number of realizations in which 
the maximum number of interferers does not exceed N (the maximum number of 
interferers that can be tolerated at any instant of time during the tagged packet trans- 
mission) and is given by [Ref. 18] 



where C,{n) is the number of realizations in which the maximum number of interferers 
equals j. 

The capture effect in a random access channel occurs naturally when fading reduces 
the joint power of interfering packets, as reported in (Ref. 19] for a slotted Aloha channel 
in a uniform Rayleigh fading environment. This thesis extends these results to unslotted 
Aloha in a Rician fading environment. 

It is assumed that a tagged packet will capture the receiver if its signal-to- 
interference ratio exceeds a threshold y 0 . 

The probability density function of a packet's power in a Rician fading channel is 




(4.2) 



where 





(4.3) 



The same expression for a normalized mean power a = 1 is 



f P (p) = exp( -Aq)I 0 (2A 0 Jp ) exp( -p) 



(4.4) 
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Thus given n interfering packets during the transmission of a tagged packet, the 
density function of the power of n interfering packets is [Ref. 17] 

n 

SrW«) = -ZT z c A n WMf (4.5) 

My] 

Substituting equation (4.4) into equation (4.5) we get 



n 

m») = 4- z c /")( ex P( -AlMUjy) “P( (4.6) 

^ : f n , 



where j® denotes the j-fold convolution. 
Appendix B illustrates the derivation of 



[ exp( - Ao)I 0 (2A 0x fy ) exp( -y)} /0 for j = 1,2, 3, 4, 5 and 6 



The density function of the power of the signal-to-interference ratio Z = 
tained by 




is ob- 



r oo 

fziAn) = yfx<yz]fy<y\n)dy 
J o 



(4.7) 



Thus 



fz( 2 I") = I T exp( - Aq)J 0 (2AqJ} ! z ) exp( -yz)\ Y C(«)[ exp - AqI 0 {2A 0 Jy ' )• 

2 J ° Uk2, U 

■ exp( — J 7 )} 70 } 



(4.8) 



The capture probability is given by 



Pr{capiure\n, /7>1) = 1 — P r {0<Z<yQ\n, w^>l) 



(4.9) 



Therefore we have 



Pr{capture\n , «>1} = 1 — 




(4.10) 



46 



Appendix B illustrates the derivation of the probabilities of capture for n = 1,2. 3, 4,5 
and 6. 

Channel throughput S is defined by 



OP 

S = GPr{capture } = G^^Pr{capture\n}Pr{n) (4.11) 

n=0 

Replacing equations (4.1) and (4.9) in equation (4.11), the throughput of unslotted 
Aloha in a Rician fading channel is given by 



S = 



oo 




> 7=0 



exp( — 2G)Pr{capture\n) 



and given that Pr{capture |0] = 1 it reduces to 



( 4 . 12 ) 



G exp( — 2G) 



OO 




;*=! 



Pr{capiure\n) 



(4.13) 



Figures 32. 33 and 34 illustrate the throughput for several capture thresholds y 0 us- 
ing A 0 = 0.1, A 0 = 0.5 and A 0 = 1 respectively. Finally, Figure 35 illustrates the through- 
put for the same capture thresholds y 0 , but letting parameter A 0 = 0 (Rayleigh fading). 
As expected the same throughput of unslotted Aloha in a Rayleigh fading environment 
[Ref. 17] were obtained. 
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Figure 32. Unslotted Aloha throughput (A 0 = 0. 1) 
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Figure 33. Unslotted Aloha throughput ( A 0 = 0.5) 



49 




Figure 34. Unslotted Aloha throughput ( A 0 = 1) 
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Figure 35. Unslotted Aloha throughput ( A 0 = 0) 
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V. CONCLUSIONS 



This thesis can be divided into two main parts. The first part describes the imple- 
mentation of a software simulator for multipath fading in a maritime mobile satellite 
channel. The second part analyzes the throughput for an unslotted Aloha channel in a 
Rician fading environment. 

The fading simulator generates Rician fading sequences for any value of frequency, 
velocity and Rice Factor (K). Fades can be simulated for communication links using a 
geostationary satellite or low altitude satellites. 

Unfortunately, it was not possible to obtain data from experiments done in several 
institutions like JPL and CRC in order to compare results and make necessary adjust- 
ments to the simulator. 

For the second part of the thesis, a set of equations are derived to compute the 
throughput of unslotted Aloha in a Rician fading channel. The obtained equations are 
tested for several values of the parameters A 0 and y 0 . Finally, a comparison is made with 
results obtained for a Rayleigh fading channel [Ref. 17] by setting the parameter A 0 equal 
to zero, and similar results are obtained. It is noted that the throughput of a Rician 
fading Aloha channel depends on the direct component power A\. As A\ increases, the 
joint interference power also increases and consequently the throughput is reduced. 
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ooooooo n on non ooooooooooo 



APPENDIX A. SOFTWARE RICIAN SIMULATOR 



This program generates a Rician fading sequence given variance, frequency and ve- 
locity. 





wV 


VARIABLE DEFINITIONS ***** 


A 


= 


LOWER LIMIT 


B 


= 


UPPER LIMIT 


FCT 


= 


FUNCTION 


S 


= 


PARTIAL SUM VALUE 


H 


= 


INTERVAL LENGTH 


AC 


= 


SIMPSON'S RULE VALUE 


N 


= 


NUMBER OF EVEN SUBINTERVAL 



NN2 , NN , I , CC = COUNTERS 

X,H,D = COMPONENTS OF THE RULES 

VARIABLE DECLARATIONS 

COMPLEX RF( 0: 1024), R(0: 1024),DSV(0: 1024) 

REAL SIC 0: 1024),SQ(0: 1024) , V,ME ,G( 0: 500) ,F( 0: 500) , Y,W,C0NST( 0: 1000 
&) ,W0 , DS 

REAL* 8 FR,WS ,WL , A , B ,FCT, S ,X,H,D,AC(21) ,PI,Z,RGF(0: 1000), AX(0: 1000) 
& , IGF( 0: 1000) ,MAG(0: 1024), TS 
INTEGER I ,NN2,N,NN,NS 
MAIN PROGRAM 

Nl: NUMBER OF SAMPLES OF THE GAUSSIAN RANDOM GENERATOR 
N 1=5 00 

N2: NUMBER OF COEFFICIENTS CALCULATED FOR THE SHAPING FILTER 
N2=21 

NP = Nl + N2 

GAUSSIAN RANDOM GENERATOR 1 WITH VARIANCE 1, MEAN 0 
ISEED=23 

WRITE(*,*) 'ENTER VARIANCE' 

RE AD(*,*) V 
ME=0. 0 

DO 10 1=0, Nl 

CALL GAUSS ( ISEED,V,ME,W) 

G(I)=W 

10 CONTINUE 

GAUSSIAN RANDOM GENERATOR 2 WITH VARIANCE 1, MEAN 0 
DO 11 L=0 ,N1 

CALL GAUSS ( ISEED, V,ME ,Y) 

F(L)=Y 

11 CONTINUE 

WRITE( 8 ,*) (G(I),I=0,N1) 

WRITE (4,*) (F( I) ,I=0,N1) 

CALCULATION OF SHAPING FILTER COEFFICIENTS 



TO CALCULATE A DEFINITE INTEGRAL FOR A FUNCTION USING SIMPSON'S 
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c 

c 



12 



C 



14 



15 



13 

20 



21 



RULE. 

NS=21 

PI=DATAN( 1D00)*4 
C=2. 997925E+08 

WRITE (*,*) ‘ENTER OPERATING FREQ’ 

READ(*,*)FR 
WO = 2 * PI * FR 

WRITE (*,*) 'ENTER SPEED IN KNOTS’ 

READ(*,*)WM 
IF (WM .EQ. 0) THEN 
DO 12 K=0 , N 1 

DS = COS(WO*K) 

RF(K) = CMPLX(G(K),F(K)) 

DSV(K) = CMPLX(DS , 0. 0) 

R(K) = DSV(K) + RF(K) 

MAG(K) = CABS(R(K) ) 

WRITE(4,*) MAG(K) 

CONTINUE 

ELSE 

WL=C/FR 

WRITE(*,*) 'WAVELENGTH = ' ,WL 
WM=2*P I * ( WM* 1852/3600) /WL 
WS=WM/2/PI 

WRITE(* ,*) 'DOPPLER FREQ = ' ,WS 

B=0. 97*WM 

WS=2*WM 

TS=2*PI/WS 

WRITE (*,*) 'ENTER (EVEN) N FOR SIMPSON INTEGRATION' 
READ(*,*)N 



A=0 

H=(B-A)/N 

D=H/3D00 

WRITE( * ,*) ' ******COMPUTING: 
DO 13 K=0 , (NS-l)/2 
S=FCT(WL,K,WS) 

DO 14 1=1 jN-1 , 2 
WL=A+I*H 

S=S+4. 0*FCT(WL,K,WS) 
CONTINUE 

DO 15 1=2, N-2, 2 
WL=A+I*H 

S=S+2. 0*FCT(WL,K,WS) 
CONTINUE 
S=S+FCT(B,K,WS) 
AC(K+1+(NS-1)/2)=S*D 
CONTINUE 

DO 20 1=1 , (NS-1) /2 
AC(I)=AC(NS+1-I) 

DO 21 J=0 , (NS-1) 

AX( J)=AC( J+l) 

CONTINUE 



PLEASE WAIT*******' 



DO 25 N=0,NS-1 
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AX(N)=AX(N)*(0. 54-0. 46*DCOS(2*PI*N/(NS-l) ) ) 

25 CONTINUE 
C 

C CONVOLUTION OF GAUSSIAN RANDOM VECTOR 1 AND SHAPING FILTER 1 
DO 26 N=N2,N1+N2 
AX(N)=0. 0D+00 
26 CONTINUE 
C 

CALL CONVOL (N1,N2,G,AX,RGF) 

C 

C CONVOLUTION OF GAUSSIAN RANDOM VECTOR 2 AND SHAPING FILTER 2 
CALL CONVOL (N1 ,N2 ,F, AX, IGF) 

C 

C APLYING DOPPLER EFFECT TO DIRECT SIGNAL 
DO 30 K = 0 ,NP 
DS = COS ((WO + WS)*K*TS) 

C 

RF(K) = CMPLX(RGF(K) , IGF(K) ) 

C ADDING RAYLEIGH SIGNAL TO A CONSTANT TO OBTAIN RICEAN 
DSV(K) = CMPLX( DS , 0. 0) 

R(K) = DSV(K) + RF(K) 

WRITE( 8 ,*) R(K) 

MAG(K) = CABS(R(K) ) 

WRITE ( 4,*) MAG(K) 

30 CONTINUE 
END IF 
STOP 
END 
C 

REAL* 8 FUNCTION FCT(WL,K,WS) 

REAL* 8 WL,WS,PI,D 
INTEGER K 

PI=DATAN( 1D00)*4D00 
D= 2 D 0 0 * P I * K* W L / W S 

FCT=2D00/WS*( IDO 0 - ( 2D00*WL/WS ) **2 ) **( -0. 25)*DCOS(D) 

RETURN 

END 
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APPENDIX B. DERIVATION OF FORMULAS FOR UNSLOTTED 

ALOHA THROUGHPUT 



A. J-FOLD CONVOLUTION 

Given the expression 

[e- Al k(2A^)e- y ~\ (B. 1 ) 

in order to compute the j-fold convolution for j = 1,2, 3, 4, 5 and 6, we evaluate the 
characteristic function and obtain 






e^UJA^e-^dy 



(B.2) 



Replacing the modified Bessel function of order 0 by its series approximation to the 
fourth power we have 



<!>>{«) = 



-AS 



e-y^y + e-'A^ + m -£-/e iuy + 4^/^ 



. — y ^ 0 4 Jux 

+ e ■ — -ry e 
0/6 v 



dy 



(5.3) 



This is equivalent to 



= e 



-4 



OO 

e 

— OO 

~oo 






dy + Al\ y 

— OO 



-0-P)yj , 

ye dy + — 



'OO 



y 2 e-°- Ju)y dy 



+ 



36 



. 8 poo 



+ ^- | Se-^»Jy 



(B.4) 



After evaluating the integrals we obtain 
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O y{u) = e 



-4 



+ A o 2 



+ 



(!-» (I -juf 4 (1 -ju) 



+ 



A 



36 n _ 



(1 -juf 



+ -rf 



24 



576 (1 -juf 



(B.5) 



Now, the j-fold convolution in time domain is equivalent to the jth power in the 
frequency domain, that is 



Evaluating the different jth powers and taking the inverse transform to each expression 
we obtain the desired j-fold convolutions. 

For j = 1 



10 “A 



(fyW^ = e 



e y + 2A\ye y -1- 0.5 — -y 2 e y -I- 0.1666 -~y 3 e y 



y. 



(5.6) 



For j = 2 

,.-n 2 ® _ „-i4 



m or = *■ 



j 7 j 0 

“ y i ^ 0 2 —y . J 0 3 —y . i ■*» -> ^ 0 4 — _y 

' + 2 — y e + 2 — >- <? ' + 1.3j 33 — y e 



3! 

■ 10 



4 S AO A2 

+ 0.5S33 4f A^ + 0.1666 -y-y 6 e~ y + 0.0277 -%-vA 



O' 
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7! - 



For j = 3 



6 



tr t \\ 3 ® -3 An 1 2 -y Aq 3 _y , „ f Aq 4 _ y Aq 5 _ v 

(w)) = « 0 —y e 7 + o — J 7 * + 4.5 —y e y + 4.0 — y e ' 



A* 

\ 7 0 6 — 

+ 3.25— y e • 



j 8 j 10- ) j 12 

^/e~ y + 1.75 A~* + 0.7083 A”" 

+ 0.2083 y 9 e- y + 0.04\6^-y ]0 e~ y + 0.0046 -yyr 



(5.7) 



(5.8) 
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For j = 4 



(/Kk )) 40 = e 4A ° ~^y 2e y + 4 ~y- y * e y + 8 ~y-y 5e y + 10,6666 ~&~y 6e ' 

A d 0 ) 



+ 9.6666 ~~y 1 e~ y + 8 



14 



10 



A 12 A 

—y e + 4 ~gry e + 2 -> 3 -> 3 "joT^ 1 

A6 AS ,20 

+ 0.8958 - { jry U e~ y + 0.2685 -y^ 1 2 



A, 

+ 0.0092 — 



22 

! 0 14 -> 



A n A 20 

0 12 — y ia/'aio 0 13 — y 

"j -yry e y + 0.6018 ~Yy-y e 



For j = 5 






,5® . 



-~ 5 A q 
, e o 



j 2 i 4 j 6 

1 4 — y . c 5 — y ioc A 0 6 — y o-.-*-. A 0 7 - 

— y * + 5 * + 12.5 —y e + 20.8j3j —y e 

a 8 j 1 0\ J 12 

+ 25.8333 -r?- vV* + 25.1666 y\~ y + 19.861 1 -^ry W e~ y 
o: y. 1U! 

j 14 j 16 ,18 

+ 12.9166^y n £^ + 6.9791 2.6157 

>2° i 22 ,24 

+ 1.1655 -y^py 1 V* + 0.353 ^-y ]: 5 e~ y + 0.0848 

A 26 a 28 ^30 

+ 0.0254 + 0.00096 ~ry l *e~ y + 0.000128 ^ry ]9 e~ y 

17! ^ 18! * 19! * 



-y 

(B.9) 



•y 



(B. 10 ) 
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For j = 6 



(/VOO)' 



,6© 



<T^o 



2 

1 5 — y . s '^0 6 — y . , Q ^0 7 — y . ^ r ^0 8 • 

—y e y +6—y e y + 18— y e + 36— y e 



0 8 -v 



A S J 10 l A 12 

c ^ h c 0 9 — y t s *s r 0 / 10 — y , /i /i i / ^ 0 11 — y 

+ 5 - ) - 75 “pi - y e + 6^-5 “ioT^ e + 61.4166 -yyr> e 

+ 49.6666 -jjry 12 ^ + 34.0208 + 19.8842 

+ 9.0069 -%-y 5 e~ y + 3 ^-ry x6 e~ y + 0.9184 ^%-y ]1 e~ y 
lo! 16! 17! 

A 26 A 28 A 30 

+ 0.3322 -jjry 18 e~ y + 0.1013 + 0.0181 

+ 0.0036 4rr v^'e^ + 0.000385 4rr A- J 



(All) 



20 -> 



B. PROBABILITY OF CAPTURE 



First, evaluating equation (4. S), f z (z \ n) for n = 1 2, 3, 4, 5 and 6, the following ex- 
pressions are obtained 



For n = 1 



fz( z I 1) 



C, (\)e- ul Y 



Iv 



1 



°=1 
o — 1 



4i*~ 

£=3 



(1+2) 

gig- l)(g- 2) 



5+1 



+ Ar 



j 



g(g~ 1) , ,4 

+ A 0 . 



o=7 

o — 



(l+2)" +1 



3 

J 6 _2 — 

0 36 



+ An -77- ^Kg - 3 
g=4 



(1+2) 



g- 



5+1 



-i (A 12) 



(g — 4)!(1 + zf +1 
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For n = 2 



/z(z I 2) 



For n = 



fzb I 3 ) 
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For n = 6 
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where the Cj(n) values are defined in equation (4.3), and K v K h , K„ K jt K m , and K u are given 
in Table 3. 
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Table 3. CONSTANT VALUES USED IN EQUATIONS (B.12) THROUGH 
(B.23) 
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Now, the probabilities of capture are evaluated by replacing equations (B.12) 
through (B.17) into equation (4.10), and we get 
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For n = 2 
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For n = 3 
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where y 0 is the preset threshold and the values of A^, A),, A!), A), A4, andA^ are given 
in table 3. 
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